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Abstract. Extreme-Mass- Ratio Inspirals (EMRIs) are one of the most promising sources 
of gravitational waves (GWs) for space-based detectors like the Laser Interferometer Space 
O Antenna (LISA). EMRIs consist of a compact stellar object orbiting around a massive black 

hole (MBH). Since EMRI signals are expected to be long lasting (containing of the order of 
j.^ ■ hundred thousand cycles), they will encode the structure of the MBH gravitational potential in 

a precise way such that features depending on the theory of gravity governing the system may 
be distinguished. That is, EMRI signals may be used to test gravity and the geometry of black 
holes. However, the development of a practical methodology for computing the generation and 
| propagation of GWs from EMRIs in theories of gravity different than General Relativity (GR) 

, has only recently begun. In this paper, we present a parameter estimation study of EMRIs in 

£S1 ' a particular modification of GR, which is described by a four-dimensional Chern-Simons (CS) 

CO , gravitational term. We focus on determining to what extent a space-based GW observatory like 

LISA could distinguish between GR and CS gravity through the detection of GWs from EMRIs. 

o 

CN , 

1. Introduction 

It is known that surrounding the central MBH of many quiescent galaxies there are around 
10 7 — 10 s stars forming a cusp/core (see e.g. p]). As a consequence of mass segregation and 
large scattering encounters between stars in the cluster, the largest stellar-mass compact objects 
(SCOs), like neutron stars, m ~ 1.4Mq, white dwarfs, m ~ 0.6Mq, and stellar-mass black holes, 
m ~ 1 — 50M Q , may be perturbed onto an orbit that passes close enough to the central MBH 
and become gravitationally bound to it. It is expected that the capture of a SCO by a MBH 
may be a frequent phenomenon in the Universe. The MBHs of interest for a GW observatory 
like LISA [2] have masses in the range M, ~ 10 5 - 1O 7 M . In addition, there are observations 
suggesting the existence of Intermediate-Mass Black Holes (IMBHs), whose masses are in the 
range M, ~ 10 2 — 10 4 M Q (see e.g. [3]). These may form EMRI-like systems, either an IMBH- 
MBH inspiral or an SCO-IMBH system. Space-based detectors like LISA will be sensitive to the 
inspiral of an IMBH into a MBH, whereas ground-based detectors like the Einstein Telescope [1] 
will be sensitive to the inspiral of SCOs into IMBHs. 

Once the SCO is captured, the binary loses energy and angular momentum through the 
emission of GWs. This drives the inspiral of the SCO until it plunges into the MBH. Due to 
the fact that the mass of the SCO is much smaller than the mass of the MBH, with typical 
mass ratios lying in the range \i = m/M, ~ 10~ 7 — 10~ 3 , the inspiral can be modelled as a slow 



adiabatic process in which the orbital time scale is much shorter than the timescale on which 
the orbit shrinks. 

During the last year of evolution before plunge, the orbit of the SCO spends of the order of 
10 5 cycles in the strong field region of the MBH, mapping its geometric structure (which can 
be characterized in terms of mass and current multipolar moments). This map is encoded in 
the shape and timing of the GWs emitted [5]. If we can decode this map, EMRI detections will 
allow us to study different aspects of the MBHs located at the centers of galaxies, the stellar 
dynamics around galactic centers, etc. In particular, EMRI GW signals may be very useful to 
test the spacetime geometry of MBHs El El El ESI EH E2 E3] and even alternative theories of 
gravity (see, e.g. [HI ESI E3 ESI E3 ES] ) , which is the subject of this article. 

Due to the myriad of alternative theories of gravity available, the question that arises here 
is which kind of theory do we choose to test? Alternatively, what kind of phenomena can we 
expect to test using EMRI observations? In this paper, we will focus on Dynamical Chern-Simons 
Modified Gravity (DCSMG) which is a 4D effective theory of gravity that arises in certain low- 
energy limits of string theory and can also arise in Loop Quantum Gravity (see |19j for a review of 
DCSMG). DCSMG has been proposed as an explanation for the cosmic baryon asymmetry |20j 
and makes predictions for the polarisation of the cosmic microwave background [21J. In addition, 
it has potentially observable signatures in EMRI GW observations |22j . 

In this paper we summarize the necessary techniques that we have developed to estimate 
the capability of a LISA-like GW detector to distinguish between GR and DCSMG using GW 
signals from EMRIs, extending and improving the results of [22J. A more detailed account of 
these developments will appear elsewhere [23] . Due to the large difference between the masses 
of the components in the binary, an EMRI signal can be modelled accurately using black hole 
perturbation theory (see e.g. [24]). in which the SCO is represented as a particle orbiting in 
a given MBH spacetime background. In this framework, at the lowest order of approximation 
the orbit of the particle is a geodesic of the MBH spacetime, but radiation reaction (RR) 
effects, i.e. the effects that arise from the interaction of the SCO with its own gravitational field, 
accelerate the trajectory and then drive the inspiral of the SCO. Depending on the theory of 
gravity that describes the system, the RR effects drive the inspiral in different ways . 

This paper is organised as follows. In Sec. [2] we summarize the formulation of our EMRI 
model in DCSMG. In Sec. [3] we review the basics of signal analysis employed in our study. 
Finally, in Sec. HI we present results concerning whether a GW detector like LISA will be able 
to discriminate between an EMRI in GR and one in DCSMG. 

Throughout this paper, we follow the conventions of Misner, Thorne and Wheeler [25] : Greek 
letters will represent spacetime indexes and Latin letters in the middle of the alphabet 
will represent spatial indexes only. Partial derivatives are denoted by d a h = dh/dx a = h j(X , 
while, for any quantity h, covariant derivatives are denoted by V a /i = h. a . In addition, we use 
the Einstein summation convention and geometrised units, i.e. G = c = 1 and the gravitational 
constant is given by k = 1 /(167rG). Finally, we use a bar, e.g. h, to denote background quantities. 



2. Formulation of our EMRI model 

We approximate our system following the procedure described in [22J, in which an EMRI 
is modelled as a binary composed of a point particle moving in a modified Kerr MBH 
background geometry. That is, the MBH geometry is slightly modified as a consequence of 
the CS gravitational modifications to GR. The motion of the particle, at the lowest order of 
approximation, is given by a geodesic of this CS-modified MBH geometry. The action describing 
the system is given by: 

a /3 

S = KSEH + -jScS + 7T<Sl9 + ^mat , (1) 



where the different terms are: 
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which represent the Einstein-Hilbert action [Eq. (|2J)], the CS gravitational correction [Eq. ©], 
the CS scalar field action [Eq. 0], and the action of the matter degrees of freedom [Eq. fl5J)]. 
The coupling constant, a, is associated with the Pontryagin density, *RR := ±e* SflL 'R% eta ^R' 3 5 , 
which is coupled to the CS scalar field, The coupling constant, /3, determines the strength of 
the coupling of the CS scalar field with the stress-energy distribution. 

Following [26j, we can construct the modified MBH geometry and the CS scalar field using the 
slow-rotation (a/M, <C 1, where M, and a = |S/M.| are the MBH mass and spin respectively) 
and weak-coupling approximation. Then, the spacetime metric, g^„, and the CS scalar field, 
can be written as the sum of a background and a perturbative quantity characterized by the 
perturbative parameter v: 
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The MBH metric has the following non-zero components in Boyer-Lindquist-like coordinates 
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where S = (r 2 + a 2 ) 2 — o 2 Asin 2 6, p 2 
The CS parameter £ is given by 
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Notice that the only metric component that contains the CS parameter £ is g t( p and that it 
behaves as r -4 for M./r <C 1, decaying faster than the rest of metric contributions and becoming 
negligible at large distances. Under the preceding approximations, the CS scalar field is given 
by 
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The numerical evolution of the system comprises two stages, namely the computation of the 
SCO trajectory around the MBH and then the computation of the corresponding waveform 
based on this trajectory. To generate the trajectory, it has been shown in [22] that in the 
test-mass limit the SCO must follow geodesies of the CS deformed MBH background. It turns 
out that the metric [Eqs. ([8])- ([12p ] has almost identical symmetry properties as the Kerr metric 
(at least at the level of approximation used for its derivation). As a consequence, the geodesic 
motion of the particle is obtained by solving the following set of Ordinary Differential Equations 
(ODEs): 

i = i K + L z 5g% s , (15) 

$ = 4>K-E8gf, (16) 

r 2 = r 2 K + 2EL z f5g$ s , (17) 

e 2 = 2 K, (18) 

where E is the energy and L z is the z-component of the orbital angular momentum of the 
particle per unit mass. Here, a dot denotes differentiation with respect to the proper time and 
(tx, ^k^Ri 4>k) are the corresponding expressions for geodesies in a Kerr MBH background: 

P H K = - (a 2 E sin 2 9 - aL z ) + (r 2 + a 2 ) ((r 2 + a 2 ) | - ^ V (19) 
P 2 4>K = -(aE--^-)+(r 2 + a^-^, (20) 
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P 4 r 2 K = [(r 2 + a 2 )E-aL z ] 2 -A[Q+(aE-L z ) 2 + r 2 ] , (21) 
p*d 2 K = Q - cot 2 9L 2 Z - a 2 cos 2 9 {1 - E 2 ), (22) 



where Q denotes the Carter constant per unit mass squared. The correction 5g^ s appearing in 
Eqs. (|15|) - ([18| ) is the CS correction, which is given by: 

5g ° S = vn?R ( 70r2 + 120rM * + 189M « ) • (23) 

In order to avoid the presence of turning points in the integration of these equations of motion, 
we introduce two new angle coordinates, ip and x> associated with the radial and polar motion 
respectively, which are defined by the equations: 

r = -x — T ' cos2 = cos2 ^in cos2 X ' ( 24 ) 
1 + e cos ip 

where p and e are the semilatus rectum and the eccentricity of the orbit respectively, and 9 min 
is the minimum value of the polar angle 8 in the orbit. This value is related to the angle of 
inclination of the orbit, 9 inc (see Table [1]), through 

9 . =--6. . (25) 

mm 2 mc V / 

A given orbit can be characterized in terms of the constants of motion (E, L z , Q) or, alternatively, 
in terms of the orbital parameters (e,p, 9 inc ). To transform between these two sets of constants 
we use the equations at the turning points of the radial and polar coordinates. At 9 = 9 min we 
have that 9 = 0, leading to an expression for Q [22 : 



Q = cos 2 6> min 
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The additional two equations needed are those that define the radial turning points, that is, the 
pericenter, r pcri , and apocenter, r apo , which we relate to (e,p) in the usual way: 
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Eq. (j26|) together with the equations r(r per i) = r{r apo ) = allow us to relate the values of the 
constants (E,L Z , Q) to those of the orbital parameters (e,p,6 inc ). 

Until now we have only considered geodesic motion. Constructing waveforms from geodesic 
trajectories can lead to a a confusion problem in the sense that it is possible that an EMRI 
waveform in GR matches an EMRI waveform in DCSMG but with different parameters, given 
that the motion in both cases is going to be characterized by just three fundamental frequencies 
of motion. This degeneracy is overcome when RR effects are taken into account, since these 
tend to break such confusion, as the RR effects drive the inspiral in different ways depending 
on the theory of gravity that describes the EMRI system. Nevertheless, it was shown in [22] 
that, at leading order, the GW emission in DCSMG takes the same form as in GR. So, as a first 
approximation we describe RR effects in DCSMG using expressions from GR. In particular we 
use the fluxes for the constants of motion (E,L Z , Q) provided in [27]. In principle one might 
think that this would lead to exactly the same inspiral evolution in DCSMG and GR, but it is 
not the case due to the different relations between the constants of motion (E,L Z , Q) and the 
orbital parameters (e,p,8 inc ). 

To compute the EMRI waveforms in DCSMG, we first solve the ODEs given in Eqs. (|15|) - 
(fTBj) for the angle variables (ip(t), x(t), ^f(t)) (see e.g. [2H])- The scheme for evolving a particular 
EMRI in DCSMG can be summarized as follows: For a given set of initial orbital parameters, 
( e (°) , p(°) , QyJ ) , we find the associated initial constants of the motion, (E^ , L z , Q(°) ) , which are 
different from the ones that one would obtain in GR for the same orbital parameters. Then, the 
time evolution of the "constants" of motion, (E^ \L Z °\ Q^), the RR effect, is computed using 
the fluxes of [27]. From the values of the constants of motion, (E^°\L z °\ Q^), their evolution 

due to RR, (E^°\ L z °\ Q^), and the value of the geodesic radial period, T r (the coordinate time 
to go from the apocenter to the pericenter and back again to apocenter), we obtain the values 

of the new constants of motion, (E^\L Z ; Q^)- Finally, from the new values of the constants 
of motion we find the new values of the orbital parameters {p^ l \ e^ l \ 6^})- This algorithm 
is iterated along the whole EMRI evolution to obtain the SCO orbit around the MBH. The 
gravitational waveform is then generated using the multipolar formalism given by Thorne [29J 
up to quadrupolar order, by identifying the Boyer-Lindquist coordinates with spherical polar 
coordinates in a pseudo flat-space. 

3. Signal analysis formalism 

In this section we briefly review the formalism of signal analysis that we use in this paper. For 
our calculations we use the LISA noise sensitivity curve (see e.g. [30]) to represent a possible 
future space-based GW detector. At low- frequencies the response of the LISA detector can be 
decomposed as the response of a pair of independent, right-angle Michelson detectors, which 
we label by a = /, // in the standard way. Each of these data streams is a time series, 
s a (t), which is a superposition of a possible gravitational wave signal h a (t) and noise n a (t), 
s a(t) = ^ a (t) + n a (t). Under the assumptions that the noise in the detector is stationary, 
Gaussian and uncorrelated between the two data streams, and that the spectral density of the 
noise is the same in the two channels, the Fourier components of the noise, n a (f), would have 
a Gaussian probability distribution and 
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where a tilde denotes the Fourier transform, an asterisk denotes complex conjugation, angle 
brackets denote the ensemble expectation value, and S n (f) is the one-sided noise spectral density 
of the data stream, which is the same for both channels. 

The vector space of signals has a natural inner product associated with S n (f): 

(a|b)=2E rmi+MM d/ (29) 

The assumption of Gaussianity then means that the probability that the noise, n, has a particular 
realisation, no, has the form 

p(n = n ) oce- (no|no)/2 , (30) 

and hence the likelihood that a signal h is present in the data stream s is proportional to 
exp(— (s — h|s — h)/2). The 'best-fit' waveform will be the one that maximises (s|h) and, thus, 
provides the maximum likelihood parameter estimate. The expected signal-to-noise ratio (SNR) 
when filtering with the correct waveform is 

(h|h) 



SNR= v, n = vW- ( 31 ) 
rms(h|n) 

where the second equality follows from the fact that the expectation value of (a|n)(b|n) is 
(a|b) [SI]. 

A given template family h(i; A) will depend on a set of parameters A. If we expand 
exp(— (s — h|s — h)/2) around the best-fit parameters, Ao, by writing A = Ao + SX, we can 
derive a probability distribution function for the parameter uncertainties 



p(SX) = M 



exp (-~T jk S\H\ k \ , (32) 



where W = v / det(r/2vr) is the normalization factor and is the Fisher information matrix 

(33) 
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The variance-covariance matrix for the waveform parameters is therefore given by the inverse of 
the Fisher matrix 

(5\i5\ k ) = (r- 1 ) jk + 0(1/SNR) (34) 

and so we can estimate the precision with which we will be able to measure a particular 
parameter, Aj, using LISA observations by computing the component T^ 1 of this inverse matrix. 

In our simulations, the space of EMRI waveforms in DCSMG, {h(t;X l )\i = l,...,N} is 
characterised by the following 15 parameters (see Table [1] for a description of these parameters): 

A = {M„a,fi,e ,p ,e^ c ,i,e s ,4> s ,e K ,(t) K ,D L ,i; ,Xo,4>o} ■ ( 35 ) 

The Fisher matrices for EMRI waveforms have very large condition numbers (the ratio of the 
largest to the smallest eigenvalues) and so we use an LU decomposition to invert them, where 
the Fisher matrix [Eq. (|34p ] is written as the product of a lower triangular matrix and an upper 
triangular matrix [33J. In addition, to verify that the error estimates that we obtain are reliable, 
we employ the maximum-mismatch criterion proposed in 



4. Results 

The main goal of our study is to determine how well a GW detector like LISA could estimate 
the parameters of an EMRI system (see Table [1]) in DCSMG and in particular how well it could 
constrain the CS parameter £ and hence distinguish between an inspiral in DCSMG and one in 
GR. To that end, we have followed the approach of [22] to model the motion of the SCO in the 
modified DCSMG MBH geometry, but this time we have included RR effects. To evolve the 
inspiral we have followed the algorithmic scheme described above. Once the trajectory of the 
SCO has been computed, we obtain the waveforms by using the multipolar formalism described 
by Thorne in but only keeping the quadrupolar term. To compute the Fisher matrix, 
the waveform derivatives, <9jh = dh/dX 1 , are obtained numerically through a five-point finite 
difference rule. Finally, the response of the detector and the Doppler modulation for a detector 
like LISA are computed using the expressions in [30J and in [35J . 



Table 1. Summary of the EMRI parameter space employed in our analysis. The angles (9c, (f>g) 
and {Ox^k) are spherical coordinates and t stands for the initial time in the computation. 
For the parameters with physical dimensions, we indicate the units used in round brackets. 



M. 




MBH mass. 


a = 


\S\/M. {Ml) 


MBH spin. 


V = 


m/M, 


Mass-ratio. 


e o 




Eccentricity of the particle orbit at t . 


Po 




Semilatus rectum at t . 


op 




Inclination of the orbit at t . 


£ = 


^ s (m. 6 ) 


Product of the Chern-Simons parameter and the spin. 






Polar angle with respect to the ecliptic. 


<i>s 




Azimuthal angle with respect to the ecliptic. 






Spin polar angle with respect to the ecliptic. 


<t>K 




Spin azimuthal angle with respect to the ecliptic. 


D L 


(Gpc) 


Distance to the EMRI from the ecliptic baricentre. 


% 




Phase angle associated with the radial motion at t Q . 


Xo 




Phase angle associated with the polar motion at t . 


00 




Boyer-Lindquist azimuthal coordinate, (j), at t Q . 



Table 2. Parameters of the different systems used in our work, where we have fixed p = 10, 
6 inc = 0.85, 6 S = 1.1, <t> s = 0.3, 6 K = 1.4, <t> K = 0.25, D L = lCpc, % = 0.25, Xo = I, and 
6 n = 0.1. 



-> <Ps = 


0.3, e K = 


1.4, <t> K 


= 0.25, 


D L 


= lGpc, % 


System M, 


a/M, 


i 1 




eo 




A 


5 x 10 5 


0.5 


2 x 10" 


-5 


0.5 


10~ 2 


B 


5 x 10 5 


0.5 


2 x 10' 


-5 


0.5 





C 


5 x 10 & 


0.25 


2 x 10" 


-5 


0.5 


10~ 2 


D 


5 x 10 & 


0.25 


2 x 10" 


-5 


0.5 





E 


10 e 


0.25 


10~ 5 




0.32 


10" 2 


F 


10 e 


0.25 


10" b 




0.32 






We have computed parameter measurement error estimates for different systems in DCSMG 
and the analogous systems in GR (that is, setting £ = = 5£). In Table [21 we show the 



parameters of the different (systems) used in this analysis, where the distance to the source is 
given in Gpc, the angles are given in radians and the rest of the parameters are given in MBH 
mass units, while the MBH mass itself is given in solar masses. The direction of the MBH spin 
is taken to be parallel to the z-direction of the EMRI reference frame. 

As a test of our numerical set-up, in Table [3] we compare our numerical estimates in GR with 
the ones obtained in [30J. Since the parameter space employed in our analysis differs with the 
one chosen in [3D], we just compare the results obtained for the parameters common to both 
works and, as we can see from Table [3J the two methods produce very similar results. However, 
as it is shown in the table, our estimate for the error in the MBH spin is two orders of magnitude 
smaller than the one obtained in [30J. In our analysis, we include the effect of the MBH spin 
directly in the geodesic equations, whereas in [30] the effect of the MBH spin is introduced as a 
post-Newtonian correction in their post-Keplerian equations of motion. Our analysis therefore 
includes higher order spin corrections which may explain some of the improvement that our 
results show. A similar improvement in spin measurement precision was observed in [33J who 
used a similar method based on exact Kerr geodesic trajectories. 

Table 3. Parameter accuracy estimates for the inspiral of a 10M Q SCO into a 10 6 M Q MBH. 
It is assumed that we have observed the last year of inspiral and that the system is an EMRI 
described by GR, i.e., £ = 0. We compare the results we have obtained with the ones obtained 
for the same parameters by Barack and Cutler |30j . 



Parameter 


GR (this work) 


GR Q30J) 


AlnM. 


8.6 x 1(T 4 


9.2 x 10 -4 


Aa/M. (AT.) 


9.5 x 1(T 6 


6.3 x 1(T 4 


Aln/i 


1.4 x 1(T 4 


9.2 x 1(T 5 


Ae 


2.7 x 1(T 5 


2.8 x 1(T 4 


AlnDx^ , Aln(»/D L ) W 


6.5 x 10 -2 ^ 


3.7 x 1(T 2 W 



We have considered two different kind of EMRIs, namely one composed of a MBH with 
M, = 5 x 1O 5 M and a SCO with mass m = 10M@, and the other one composed of a MBH 
with mass M, = 10 6 Mq and a SCO with mass m = 1OM . For the former, we consider an 
observation lasting for the six months before plunge (SNR= 35) and for the latter we consider an 
observation lasting one year before plunge (SNR= 40). In Table 01 we show the error estimates 
obtained in our simulations. We typically find that an EMRI observation can determine the 
extrinsic parameters, namely A ex trinsic = {#g> <t>s> ®ki •Afo D L , tp , Xoi 0ol > t° a fractional error 
of ~ 10~ 2 , while the intrinsic parameters: Ai ntr i ns i c = {M«, a, /i, e , 0q c }, are typically measured 
with a fractional error of ~ 10~ 4 , 10 -5 , 10 -5 , 10 -5 , 10 -5 respectively. For the CS parameter £, 
we find that in the best case it can be measured with a fractional error ~ 10 _1 . These results 
suggest that LISA might be able to distinguish between GR and DCSMG using certain EMRI 
observations. Notice that in this paper we are using the concept of the LISA mission as a 
reference and, as it is shown in [23], the results obtained are also valid for the updated eLISA or 
NGO (New Gravitational- Waves Observatory) [36] by normalising to a fixed SNR. In addition, 
we find that for a given value of £, it will be better estimated for systems with smaller MBH 
mass, M m , and spin, a. This is not surprising, since the leading-order term that contains the CS 
parameter in the DCSMG metric [Eq. (|12p ] is ~ i/M®(M./r) 6 . By decreasing the MBH mass 
with fixed £, this term becomes larger at fixed r/M, and by decreasing the MBH spin this term 
is larger relative to the Kerr g t( k metric coefficient. 

Due to the relatively high computational cost of these calculations, a full (i.e., including all the 



Table 4. Error estimates for half a year of evolution (SNR= 35) for Systems A, B, C, and D; 
and for one year of evolution (SNR= 40) for Systems E and F. 



Parameter A B C D E F 



AlnM. 


1. 


.0 


X 


1(T 


-4 


9. 


.8 


X 


io- 


-3 


7. 


.8 


X 


io- 


-l 


7. 


.6 


X 


io- 


-4 


8. 


.6 


X 


io- 


-1 


8. 


.4 


X 


io- 


-4 


Aa (M.) 


1. 


.5 


X 


10" 


_5 


9. 


.5 


X 


io- 


-6 


1. 


.5 


X 


io- 


_(-; 


6. 


.3 


X 


io- 


-6 


1. 


.8 


X 


io- 


_5 


9. 


1 


X 


io- 


_(-; 


Aln^i 


1. 


.8 


X 


io- 


-4 


3. 


9 


X 


io- 


-5 


5. 


.7 


X 


io- 


-5 


2. 


.4 


X 


io- 


-5 


3. 


.2 


X 


io- 


-4 


1. 


.2 


X 


io- 


-4 


Ae 


9. 


2 


X 


io- 


-5 


9. 


2 


X 


io- 


-5 


8. 


.4 


X 


io- 


-5 


8. 


.4 


X 


io- 


-5 


1. 


.9 


X 


io- 


-5 


1. 


.7 


X 


io- 


-5 


A9£ c 


2. 


.8 


X 


io- 


-5 


1. 


.7 


X 


io- 


-5 


2. 


.6 


X 


io- 


-5 


2. 


.2 


X 


io- 


-5 


8. 


.1 


X 


io- 


-5 


3. 


.2 


X 


io- 


-5 


Am£ 




2.1 
















4. 


.4 


X 


io- 


-1 
















5. 


1 














A6 S 
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parameters) error estimation takes approximately three days using a single processor machine. 
For this reason, we have only studied a small number of points in the parameter space of all 
possible EMRI systems that we might detect. In this sense, these results could be improved by 
performing a more exhaustive study of the parameter space and also using a full Monte Carlo 
sampling. In addition, we could have considered longer gravitational waveforms, that cover the 
entire length of the LISA mission (~ 2-5 years), which would be expected to improve our ability 
to measure the system parameters. 

5. Conclusions 

In this work we have presented a first study aimed at assessing the ability of a GW detector like 
LISA to discriminate between GR and DCSMG by performing parameter estimation studies of 
EMRIs in DCSMG. We have computed the waveforms generated by a SCO orbiting in a MBH 
geometry modified by CS corrections. The parameter error estimates have been computed using 
a Fisher matrix analysis. To validate our computational framework, we have studied a typical 
EMRI system in GR and we have found agreement between our results and previous results 
found in the literature. We have applied these parameter estimation studies to estimate the 
ability of LISA to distinguish between GR and DCSMG by estimating the CS parameter £. To 
that end, we have performed simulations of the GWs emitted during the last six months before 
plunge of an EMRI system for which the gravitational wave emission is at frequencies in the 
most sensitive part of the LISA band. Our results indicate that, for certain EMRI systems, a 
detector like LISA may be able to discriminate between GR and DCSMG. We have also seen 
that the error in estimating £ decreases with the MBH mass and spin. We are now performing a 
more exhaustive study of the parameter space of possible EMRIs to better understand if these 
results hold in general [23J. 

Another subject that deserves further study is the question of the systematic errors that would 
arise by using GR waveform templates to detect EMRIs that are actually described by DCSMG. 
This can be done by using the formalism developed recently by Cutler and Vallisneri [37], which 
allows for the estimation of the magnitude of the model errors. Finally, the present work can 



be extended to other GW detectors, for instance, to study Intermediate-Mass-Ratio Inspiral 
detection with the future Einstein Telescope. 
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